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ABSTRACT 

We present the morphological analysis based on HST-NICMOS observations in the F160W 
filter (A ^ 1.6 jum) of a sample of 32 early-type galaxies (ETGs) at 1 < z < 2 with spectro- 
scopic confirmation of their redshift and spectral type. The 32 ETGs at (z) ~ 1.5 are placed 
on the Re] plane according to a Kormendy relation with the same slope of the local one 
but with a different zero-point which accounts for the evolution they undergo from z ~ 1.5-2 
to z = 0. The best fitting of their SED shows that ETGs at 1 < z < 2 are composed of two 
distinct populations, an older population (oETGs) and a younger population (yETGs) whose 
mean ages differ by about 1.5-2 Gyr. Young ETGs are not denser than local ones since they 
follow the size-mass relation of local ETGs and luminosity evolution brings them onto the 
local Kormendy and size-luminosity relations without the need of size evolution. Old ETGs 
do not follow the size-mass relation of local ETGs and luminosity evolution does not account 
for the discrepancy they show with respect to the local size-luminosity and Kormendy rela- 
tions. An increase of their R^, by a factor 2.5-3 (a density decrease by a factor 15-30) from 
z ~ 1.5 - 2 to z ~ is required to bring these galaxies onto the local scaling relations. The 
different properties and the different behaviour shown by the two populations with respect to 
the scaling relations imply different formation and evolution scenarios. The older population 
of ETGs must have formed at higher-z in a sort of dissipative gas-rich collapse ables to pro- 
duce remnants which at z ~ 2 are old and compact, scenario which can be fitted qualitatively 
by some recent hydrodynamic simulations of gas-rich mergers. Given the typical time scale 
of merging and the old age of their stellar population, oETGs should exist as they are up to 
z > 3 - 3.5. The size evolution they must experience from z ~ 2 to z ~ must leave unchanged 
their mass to not exceed the local number of high-mass (A1« > 5 ■ 10" Mq) ETGs. Thus, ma- 
jor merging cannot fit this requirement. Satellite merging, close encounters and interactions 
can help at least qualitatively in solving this problem. The younger population of ETGs can 
be formed later through subsequent episodes of merging which increased progressively their 
size and assembled their mass down to z ~ 2. At z < 2 they evolve purely in luminosity since 
episodes of major merging would bring them far from the local scaling relations. 
paolo.saracco@brera.inaf.it 
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1 INTRODUCTION 

The formation and the evolution of early-type galaxies (ETGs, el- 
liptical and bulge-dominated galaxies) occupy an important posi- 
tion among the challenges of the observational cosmology. At least 
~ 70% of the stellar mass in the local universe is locked into ETGs. 
For this reason, the understanding of their build-up and growth is 
fundamental to trace the galaxy mass assembly in the Universe. 
Their homogeneous properties, e.g. colors and scaling relations, 
make them an excellent probe to investigate the history of the stel- 
lar mass assembly of galaxies over cosmic times (see e.g. Renzini 
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et al. 2006 and references therein). Most of the recent studies based 
on samples of ETGs at z < 1 agree with considering completed 
their build up at z ~ 0.8. This statement is supported by the results 
on the evolution of the stellar mass function of galaxies which do 
not show any deficit of high-mass galaxies up to z = 0.8 - 1 (e.g. 
Fontana et al. 2004; Pozzetti et al. 2007), by the observed evolution 
of the bright end of the luminosity function of galaxies consistent 
with the pure luminosity evolution expected for early-types (e.g. 
Drory et al. 2005; Saracco et al. 2006; Zucca et al. 2006; Caputi et 
al. 2007; Cirasuolo et al. 2007; 2008), by the observed number den- 
sity of ETGs at z < 1 consistent with the one at z = (e.g. Saracco 
et al. 2005; Cimatti et al. 2006; Conselice et al. 2007) and by the 
evolution of the size-mass and size-luminosity relations compati- 
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ble with a passive luminosity evolution (e.g. Mcintosh et al. 2005). 
If on one hand the agreement among these results provide a clear 
view of the status and of the evolution of ETGs at z < I, they 
do not add stringent constraints on the mechanism(s) with which 
ETGs assemble their mass, pushing at higher z the redshift range of 
interest. 

Several studies, indeed, suggest that both the observational 
and the theoretical efforts aimed at constraining the build up of the 
stellar mass of ETGs and the shaping of their morphology should 
be focused at I < z < 2, the redshift range for which the strongest 
evolution is expected (e.g. Glazebrook et al. 2004; Amouts et al. 
2007). The picture at this redshift is far from being as clear and 
homogeneous as at z < 1 because of the difficulties in catching 
ETGs at high redshift. Indeed, to study the population of ETGs at 
1 < z < 2, a preliminary but not so obvious step has to be overcome: 
the identification of suitable samples of ETGs with secure redshift 
determination, spectral classification and confirmed morphological 
signatures similar to those of the local ETGs. In fact, up to now, 
only few samples of spectroscopically identified ETGs at z > 1 
have been collected, with no more than a tenth of galaxies mor- 
phologically confirmed each: the sample of McCarthy et al. (2004) 
resulting from the GDDS (Abraham et al. 2004) contains 10 galax- 
ies at z ~ 1.6; the sample of Longhetti et al. (2005) derived from the 
TESIS (Saracco et al. 2003) contains 10 galaxies at z ~ 1.4 and the 
sample of Cimatti et al. (2008) resulting from the K20 (Cimatti et 
al. 2002) and the GMASS surveys contains 13 galaxies at z ~ 1.6. 
On the basis of the analysis of these few samples, it is well ascer- 
tained that ETGs at z ~ 1-2 contain stellar populations formed at 
z > 3 in an intense and short lived starburst (Longhetti et al. 2005; 
Kong et al. 2006; Farrah et al. 2006; Cimatti et al. 2004; McGrath 
et al. 2007). This information provides a strong constraint on the 
star formation history of the stellar population of ETGs, that is the 
epoch at which the stars they host formed. On the other hand, this 
does not constrain how ETGs grew up: what are the time scale and 
the mechanism(s) characterizing the growing and the shaping of 
ETGs ? 

Recently, evidence for higher compactness of the ETGs at 
z > 1 with respect to the local ones have come out. Daddi et al. 
(2005) show that a large fraction of their ETGs have smaller sizes 
(efifective radii ~ 1 kpc) than local ETGs of comparable stellar 
mass, possibly implying higher stellar densities even if the presence 
of AGN in some of them could justify the compactness. However, 
other studies find similar results confirming the apparent smaller 
sizes of high-z ETGs if compared to the local ETGs with compara- 
ble stellar mass (e.g. Di Serego Alighieri et al. 2005; Trujillo et al. 
2006a; Cassata et al. 2005). These results are based on HST opti- 
cal observations sampling the blue and UV rest-frame emission of 
the galaxies and/or on seeing limited ground-based observations, 
characteristics which could affect the estimate of the effective ra- 
dius of high-z early-type galaxies. Moreover, the above results have 
been obtained by comparing galaxies at different redshift having, in 
principle, the same stellar mass. However, the stellar mass estimate 
depends on the spectrophotometric models used to fit the data and 
on the different model parameters. More recently, Longhetti et al. 
(2007) have studied the Kormendy relation for a sample of ETGs 
at z ~ 1.5 using HST-NICMOS observations. They show that these 
ETGs are at least 2 times more compact than those in the local Uni- 
verse showing that this apparent high compactness is real and not 
dependent on the wavelength of observation (see also McGrath et 
al. 2008 and Buitrago et al. 2008 for similar recent results). Some 
of these works are based on the small samples quoted above made 
of a tenth of galaxies spanning a narrow range in luminosity and 



stellar mass or on samples of candidates ETGs with no confirma- 
tion of their redshift and spectral type. 

In an attempt to provide new and stronger constraints on the 
formation of compact/dens high-z ETGs we have studied the main 
scaling relations (the Kormendy, the size-luminosity and the size- 
mass relations) for a new sample of 32 ETGs at 1 < z < 2 with 
spectroscopic confirmation of their redshift and spectral type. The 
morphological analysis for the whole sample is based on HST- 
NICMOS imaging in the F160W filter {X~ 1.6 jim) which samples 
the rest-frame R-band at the redshift of the galaxies. The sample 
spans 3 magnitudes in absolute magnitude and more than two or- 
ders of magnitude in stellar mass. This paper presents the analysis 
and the results we obtained from the study of the scaling relations 
and it is organized as follows. Section 2 is a presentation of our 
sample describing the criteria used to construct the sample and the 
data we have at hand. Section 3 describes the methodology used 
to determine the main physical properties (eff'ective radius, sur- 
face brightness, absolute magnitude, stellar mass and age) of the 
32 ETGs. Section 4 presents the Kormendy relation while Section 
5 presents a discussion of the results obtained in Sec. 4. Section 6 
shows the size-luminosity and the size-mass relations while Sec- 
tion 7 places the results in the context of the galaxy formation and 
evolution scenarios. Section 8 is a summary. 

Throughout this paper we use a standard cosmology with Hq = 
70 Km s ' Mpc"', £1„ = 0.3 and Ha = 0.7. All the magnitudes are 
in the Vega system, unless otherwise specified. 



2 SAMPLE SELECTION AND HST-NICMOS IMAGING 

The sample of early-type galaxies we constructed is composed of 
32 galaxies at 1 < z < 2 selected from different samples and sur- 
veys on the basis of their spectroscopic and morphological classi- 
fication. We restricted our selection to those galaxies having both 

deep HST-NICMOS observations in the F160W filter (A ~ 1.6 
fim) sampling the rest-frame continuum 0.55 < < 0.85 yum at 

1 < z < 2 and ii) spectroscopic confirmation of their redshift and 
spectral type. On the basis of these criteria we were able to col- 
lect a sample of 32 ETGs. The samples from which they have been 
extracted are the following: 

- 10 ETGs at 1.4 < z < 1.9 have been selected from the 
Galaxy Deep-Deep Survey sample (GDDS, Abraham et al. 2004; 
McCarthy et al. 2004). According to the spectral classification 
described in Abraham et al., we selected those galaxies having 
Class=001 (8 galaxies), i.e. pure signatures of an evolved stellar 
population, and two galaxies having signatures of a young popu- 
lation superimposed to the older one (Class=101 and Oil). These 
galaxies are listed in Tab. 1 as SA# and have been recently studied 
also by Damjanov et al. (2008). 

- 6 ETGs at z - 1.27 have been selected from the sample of Stan- 
ford et al. (1997) and belong to the cluster RDCS 0848+4453 in 
the Linx field. Their spectra show absorption features (Ca II K+K, 
Mg I and Mg II) and spectral break (B2900, D4000) similar to the 
present epoch ellipticals (Stanford et al. 2007; Van Dokkum et al. 
2003). These galaxies, listed in Tab. 1 as CIG#, have been previ- 
ously studied also by Moriondo et al. (2000). 

- 3 ETGs at 1 < z < 1.8 have been selected from the sample of 
Stanford et al. (2004) in the Hubble Deep Field-North (HDF-N) 
according to their spectral type ST< 0. 1 characterizing an old and 
passive population of stars and on their morphology. They are listed 
as HDF_#. 

- 2 ETGs (HUDF.#) at z ^ 1 .4 and at z ^ 1 .9 respectively have been 
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selected from the sample of Cimatti et al. (2008; see also Daddi et 
al. 2005) in the Hubble Ultra Deep Field (HUDF). They are clas- 
sified as early-types both on the basis of their morphology and on 
their spectral features. 

- the ETG 53W091 has been taken from the paper of Dunlop et al. 
(1996; see also Spinrad et al. 1996). The spectrum of this galaxy is 
characterized by absorption features typical of an old stellar pop- 
ulation and its light profile is bulge dominated (Waddington et al. 
2002). 

- The remaining 10 ETGs (S2F# and S7F#) come from our own 
sample of ETGs spectroscopically classified at 1.2 < z < 1.7 in the 
framework of the TESIS project (Saracco et al. 2003, 2005). The 
study of their spectro-photometric properties and of their morphol- 
ogy based on multiwavelength data and HST-NICMOS observa- 
tions are described in previous works (Longhetti et al. 2005, 2007). 

The whole sample of 32 ETGs is listed in Table 1 where we 
also report for each galaxy the photometry in different bands. All 
the magnitudes, with the exception of the F160W-band magnitude, 
are taken from the literature as quoted in the Column 4 of the table. 
The magnitude in the F160W filter is the SExtractor MAG_BEST 
magnitude (Bertin and Amouts 1997) that we estimated from the 
HST-NICMOS images we retrieved from the HST archive. All the 
magnitudes are in the Vega system. Given the different samples the 
galaxies have been extracted from, the wavelength coverage is not 
the same for all the galaxies as well as the filters used. We did not 
convert the magnitudes derived in slightly different filters from the 
different surveys to a common filter system since we would have 
introduced large and possibly systematic errors. We have preferred 
to keep the original magnitudes as given by the authors and to use 
the appropriate set of response filter functions in our analysis. We 
describe the multiwavelength coverage and the filters used in the 
various surveys in Appendix A. 

The median redshift of the sample thus collected is Zmed = 
1.45. HST-NICMOS images with the NIC2 (0.075 arcsec/pix) cam- 
era are available for the 10 galaxies of our sample and for the galaxy 
53W091, i.e. for ~ 30% of the sample. For the remaining galaxies 
the available images are based on NIC3 (0.2 arcsec/pix) camera. 
The NIC3 images relevant to the two galaxies in the HUDF were 
drizzled to 0.09 arcsec/pixel. The 1-cr limiting surface brightness 
jiiii,, of the different NICMOS images is reported in Tab. 1. The 
NICMOS mosaics for the whole sample of ETGs can be retrieved at 
the web page http : //www . brera . inaf . it/utenti/saracco/. 



3 PHYSICAL PARAMETERS OF EARLY-TYPE 
GALAXIES 

3.1 Morphological parameters 

We derived the effective radius r^ [arcsec] and the mean surface 
brightness (SB) {fi)e [mag/arcsec^] within r^ of our galaxies from 
the NICMOS images by fitting a Sersic profile (Sersic 1968) to 
the observed light profiles. The analytic expression of the adopted 
profiles is 

I(r) = I,exp{-b„[(r/r,)"" -I]] (1) 

where n = 4 and n = I values define the de Vaucouleurs (de Vau- 
couleurs 1948) and the exponential (disk) profiles respectively. We 
used Galf it software (v. 2.0.3; Peng et al. 2002) to perform the fit- 
ting to the observed profiles. The bi-dimensional Sersic model has 
been convolved with the Point Spread Function (PSF) of the NIC2 



and the NIC3 cameras modeled by means of the Tiny TinQ (v. 6.3) 
software package (Kiist 1995; Krist J. & Hook R. 2004). The fitting 
provided us with the semi-major axis of the projected elliptical 
isophote containing half of the total light and with the axial ratio 
b/a. We thus derived the circularized effective radius = -s/FJa. 

To assess the robustness and the accuracy of our estimate of 
the effective radius of galaxies, we applied the same fitting proce- 
dure to a set of simulated galaxies inserted in the real background. 
The simulations follow those described in details in Longhetti et 
al. (2007) for the NIC2 images. Here, we summarize the main fea- 
tures of the procedure followed to obtain the simulated observa- 
tions. We generated with Galf it a set of 100 galaxies described by 
a de Vaucouleurs profile with axial ratio b/a and position angle PA 
randomly assigned in the ranges 0.4 < b/a < 1 and 0<PA<180 deg 
respectively. Magnitudes in the F160W filter were assigned ran- 
domly in the range 18 < F160W < 21. Effective radii r, ,,, were 
assigned randomly in the ranges 0.2 < r^j,, < 0.5 arcsec (corre- 
sponding to 1.5-4 kpc at z ^ 1) for the NIC2 images. In order to 
verify the absence of a bias against the detection of small effective 
radii (?■(,_,„ < 0.2 arcsec, < 1.5 Kpc) in the fitting of galaxies in 
the NIC3 images characterized by a pixel scale of 0.2 arcsec/pix, 
we simulated here also a set of galaxies with 0.05 < r,_,„ < 0.2 
arcsec. The simulated galaxies have been convolved with the NIC- 
MOS PSFs and then embedded in the real NIC2 and NIC3 images. 
We used the simulations described in Longhetti et al. (2007) for the 
NIC2 galaxies and the images relevant to the sample of McCarthy 
et al. (2005) to simulate the NIC3 galaxies. We then fit the simu- 
lated galaxies with the Sersic profile and studied the behaviour of 
the resulting r^jj, checking our ability in recovering the input value 
rj ,„. In Figure 1 we plot the values of the effective radius r^ y,, ver- 
sus r^ j,, in the case of NIC2 and NIC3 images. It can be seen that 
in the NIC2 images the effective radius of the galaxies is slightly 
underestimated by Ar^"-^ ^ 0.07 arcsec on average, while in the 
NIC3 images the mean underestimate is Ar^'^^ ^ 0.03 arcsec. The 
rms measured is 0.04 arcsec in the NIC3 images and 0.02 arcsec in 
the NIC2 images and are much larger than the formal fitting error 
on the effective radius. Thus our fitting tends to slightly underesti- 
mate the effective radius r^. Fig. 1 shows also that the pixel scale of 
NIC3 does not represent a limit in the detection of small effective 
radii and thus that no bias is present in our analysis against very 
small galaxies. The small underestimate of has been taken into 
account in the following analysis and in the derivation of the mean 
surface brightnesses by adding the mean offset Ar^ quoted above 
to the best fitting values r^. In Tab. 2 we report the morphological 
parameters [arcsec] and R^ [Kpc] derived from the fitting to the 
profile of our galaxies and the mean SB in the FI60W band 

(^^yFmw ^ FmW,,, + 5log(r,) + 2.5log{2n) (2) 

where FI60W,o, is the total magnitude in this filter derived by 
galf it and reported in Tab. 2. 

The morphology for some of them, namely the CIG# galax- 
ies, the HDF# galaxies and 53W091, had been already derived in 
the rest-frame R-band from Moriondo et al. (2000), Stanford et 
al. (2004) and Waddington et al (2002) respectively. Our new esti- 
mates are in good agreement with their estimates in spite of the dif- 
ferent method used to fit the profiles. For four galaxies of our sam- 
ple, namely HUDF.472, HUDF.996, HDF.1031 and HDF.1523, 
the estimate of the effective radius derived from HST images in the 
optical bands (F814W and F850W filters) sampling the rest-frame 
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Table 1: The sample of early-type galaxies. The magnitudes are in the Vega system. 



ID 






RA 






Dec 


ref. 


NIC 


l^lira 


B 


V 


R 


I 


z 


J 




H 


F160W 


K 




S2F5.109 


03 


06 


21 


08 


-00 


17 


58 





[1] 


2 


22 


1 


24.2+0.3 


23.5+0.2 


21.8+0.1 


20.1+0.1 


- 


18.2+0 


1 


- 


17 


75+0 


05 


16.6+0 


1 


S7F5_254 


13 


34 


59 


64 


+16 


49 


10 


7 


[1] 


2 


22 


8 


25.0+9.0 


24.5+9.0 


24.0+9.0 


23.1+0.7 


- 


19.8+0 


1 


- 


19 


58+0 


09 


17.8+0 


1 


S2F057 


03 


06 


52 


25 


-00 


00 


09 


9 


[1] 


2 


22 


8 


25.0+9.0 


24.5+9.0 


23.8+0.4 


21.5+0.2 


- 


19.5+0 


1 


- 


19 


04+0 


07 


17.8+0 


1 


S2F1_389 


03 


06 


28 


03 


+00 


00 


31 


6 


[1] 


2 


22 


8 


24.3+0.5 


23.7+0.5 


23.7+0.5 


23.0+0.5 


- 


20.3+0 


2 


- 


20 


08+0 


11 


18.2+0 


1 


S2F1.511 


03 


06 


34 


04 


+00 


02 


30 


9 


[1] 


2 


22 


8 


25.0+9.0 


24.5+9.0 


24.0+9.0 


21.6+0.6 


- 


19.8+0 


1 


- 


19 


37+0 


08 


18.1+0 


1 


S2F1.142 


03 


06 


36 


51 


-00 


03 


01 





[1] 


2 


22 


8 


25.0+9.0 


24.5+9.0 


23.8+0.3 


21.5+0.2 


- 


19.6+0 


1 


- 


19 


00+0 


06 


17.8+0 


1 


S7F5.45 


13 


34 


24 


98 


+16 


45 


48 


6 


[1] 


2 


22 


8 


- 


24.2+0.4 


23.5+0.3 


22.2+0.3 


- 


19.6+0 


1 


- 


18 


59+0 


06 


17.6+0 


1 


S2F1.633 


03 


06 


35 


10 


+00 


04 


43 


6 


[1] 


2 


22 


8 


25.0+9.0 


24.5+9.0 


24.0+9.0 


22.5+0.5 


- 


20.0+0 


1 


- 


19 


32+0 


08 


18.2+0 


1 


S2F1.443 


03 


06 


31 


76 


+00 


01 


13 


4 


[1] 


2 


22 


8 


25.0+9.0 


24.5+9.0 


24.0+9.0 


23.2+0.6 


- 


20.5+0 


1 


- 


19 


69+0 


10 


18.4+0 


1 


S2F1_527 


03 


06 


43 


34 


+00 


02 


44 


8 


[1] 


2 


22 


8 


25.0+9.0 


24.5+9.0 


24.0+9.0 


22.6+0.4 


- 


20.4+0 


2 


- 


19 


80+0 


10 


18.3+0 


1 


SA12-5592 


12 


05 


22 


13 


-07 


24 


32 


6 


[2] 


3 


24 


4 


26.79+0.4 


25.00+0.09 


24.11+0.04 


23.24+0.04 


23.74+0.2 


- 




20.01+0.2 


20 


40+0 


16 


19.42+0 


20 


SA12-5869 


12 


05 


21 


55 


-07 


24 


09 


4 


[2] 


3 


24 


4 


26.90+9.0 


25.85+0.14 


24.39+0.05 


23.25+0.04 


23.78+0.2 


- 




19.25+0.1 


19 


64+0 


13 


18.58+0 


13 


SA12-6072 


12 


05 


12 


58 


-07 


23 


56 


5 


[2] 


3 


24 


4 


26.17+0.2 


25.37+0.09 


25.26+0.09 


24.10+0.08 


24.00+9.0 


- 




20.83+0.3 


21 


06+0 


19 


19.79+0 


24 


SA12-8025 


12 


05 


25 


40 


-07 


21 


24 


5 


[2] 


3 


24 


4 


26.54+0.3 


25.35+0.09 


24.55+0.05 


23.16+0.04 


22.98+0.1 


- 




19.83+0.2 


19 


94+0 


14 


18.91+0 


15 


SA12-8895 


12 


05 


14 


33 


-07 


20 


14 


5 


[2] 


3 


24 


4 


26.02+0.2 


24.67+0.05 


23.92+0.03 


23.12+0.04 


22.95+0.1 


- 




19.50+0.1 


19 


44+0 


09 


18.48+0 


12 


SA15-4367 


1-5 


23 


42 


63 


-00 


07 


11 


7 


[2] 


3 


24 


4 


26.50+9.0 


25.59+0.12 


25.03+0.07 


23.58+0.05 


22.92+0.1 


- 




21.10+9.0 


20 


75+0 


16 


19.47+0 


12 


SA15-5005 


15 


23 


45 


96 


-00 


06 


21 


3 


[2] 


3 


24 


4 


26.50+9.0 


25.09+0.09 


24.55+0.06 


23.59+0.05 


23.13+0.2 


- 




19.90+0.3 


20 


59+0 


15 


19.59+0 


14 


SA15-7543 


15 


23 


44 


83 


-00 


03 


37 


6 


[2] 


3 


24 


4 


26.50+9.0 


25.77+0.14 


24.97+0.06 


23.62+0.06 


22.73+0.1 






19.75+0.2 


19 


70+0 


10 


19.03+0 


10 


SA22-0189 


22 


17 


47 


59 


+00 


13 


27 


1 


[2] 


3 


24 


4 


99.99+9.0 


25.68+0.13 


99.99+9.99 


22.80+0.04 


22.24+0.1 






19.19+0.05 


19 


27+0 


12 


18.05+0 


11 


SA22-1983 


22 


17 


48 


41 


+00 


16 


08 


8 


[2] 


3 


24 


4 


99.99+9.0 


25.89+0.13 


99.99+9.99 


23.68+0.09 


23.28+0.2 






19.98+0.1 


20 


02+0 


19 


19.06+0 


19 


CIGJ37 


08 


48 


30 


79 


+44 


53 


34 


8 


[3,4] 


3 


24 


7 


26.3+0.1 




24.7+0.10 






21.59+0 


12 


20.4+0.09 


20 


38+0 


09 


19.53+0 


06 


CIG.65 


08 
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04 
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1 



Notes to Table 1. 



Col. 4: references for the single objects. [1] Longhetti et al. 2007; [2] Abraham et al. 2004; [3] Stanford et al. 1997; [4] Moriondo et al. 2000; [5] Stanford et al. 2004; [6] Cimatti 

et al. 2008; [10] Daddi et al. 2005; [11] Dunlop et al. 1996. Col. 5: NICMOS camera used: 2 for NIC2, 3 for NIC3. The plate scale of the NIC2 camera is 0.075 "/pixel; the 
plate scale of the NIC3 camera is 0.2 arcsec/plx. HUDF images have been drizzled to 0.09 "/pixel. Col. 6: limiting surface brightness (mag arcsec"^) at 1 — tr in the F160W 
band measured on the NICMOS images. 
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Figure 1. Comparison between the effective radius of tlie simulated galaxies 
(Tejn) and the effective radius r^y,, obtained through the fitting with the 
Sersic profile for the NIC2 (left panel) and the NIC3 (right panel) images. 



wavelength /i„s, < 4000 A, is also available from the literature. In 
the case of the two galaxies in the HDF at z ~ 1 the estimates of 
Rf we derived from the NICMOS images agree well with those de- 
rived by Van Dokkum et al. (2003) from observations in the F814W 
filter. Our estimate of for the galaxy HUDF_996 agrees with the 
one of Daddi et al. (2005) made on the F850W image while it is 
a factor two the estimate of Cimatti et al. (2008). For the galaxy 
HUDF_472 at z ~ 1.9 our estimate is a factor two both the es- 
timate of Daddi et al. (2005) and that of Cimatti et al. (2008). It 
should be noted that, given the redshift of this latter galaxy, the fil- 
ters F850W and F160W sample the profile at ^ 2900 A and 
at Arest - 5700 A respectively, two wavelength ranges which differ 
substantially for the contribution from the young and the old stars 
respectively. Thus, in this case the different estimates could reflect 
a real difference of the galaxy profile if observed at such different 
wavelengths (e.g. McGrath et al. 2008). However, we note that the 
bi-dimensional fit to the observed images of the two HUDF galax- 
ies (together with those of the two RDCS galaxies with R,, > 8 
Kpc) present significant residuals and thus their derived R, could 
be affected by large errors. The ten SA# galaxies selected from the 
GDDS have been recently studied by Damjanov et al. (2008). We 
notice that the effective radii they derive are systematically smaller 
than ours and that for 3 galaxies they derive effective radii as small 
as 0.3, 0.4 and 0.7 Kpc respectively, while we never obtain effec- 
tive radii smaller than 1 Kpc. In Fig. 2 the comparison between the 
original estimate of the effective radii R"''""" [Kpc] obtained by the 
other groups and our new estimate is shown. The different symbols 
refer to the different data sets as detailed in the caption of the figure. 
The original estimate of the effective radii [arcsec] as derived 
by the other groups are also reported in Tab. 2. 



Figure 2. Comparison between the eft'ective radius [Kpc] we obtained 
from the analysis of the NICMOS images and the one obtained by the other 
groups R"'''"^. Cyan squares are the SA# ETGs compared with the esti- 
mates of Damjanov et al. (2008); blue pentagons are the CIG# ETGs in the 
RDCS 0848+4453 compared with the estimates of Moriondo et al. (2000); 
green hexagon are the HDF# ETGs compared with the estimates of Stan- 
ford et al. (2004); red triangles are the HUDF# ETGs compared with the 
ACS-based estimates of Cimatti et al. (2008); red starred symbol is 53W091 
compared with the estimate of Weddington et al. (2002) and the black points 
are 6 out of the 10 galaxies studied in Longhetti et al. (2007) compared with 
the recent estimates of Damjanov et al. (2008). 



3.2 Absolute magnitudes, stellar masses and age 

For each galaxy of the sample we derived the R-band absolute 
magnitude Mr, the stellar mass M, and the mean age Age of the 
stellar population. We made use of the stellar population synthe- 
sis models of Chariot and Bruzual (2008, CB08) and of the best- 
fitting code hyperz (Bolzonella et al. 2000) to find the best-fitting 
template to the observed spectral energy distribution (SED) at the 
redshift of each galaxy. The set of templates considered includes 
three star formation histories (SFHs) described by an exponen- 
tially declining star formation rate 5 FR oc e"'''" with e-folding time 
r = [0.1,0.3,0.6] Gyr and two metalhcity 0.4Zo and Zq. We as- 
sumed Chabrier initial mass function (IMF; Chabrier et al. 2003). 
Extinction Ay has been considered and treated as a free parame- 
ter in the fitting. We adopted the extinction curve of Calzetti et al. 
(2000) and we allowed A^ to vary in the range ^Av < 0.6. For 24 
out of the 32 galaxies the best-fitting template is defined by SFHs 
with T < 0.3 Gyr and Ay < 0.3. 

The R-band absolute magnitude has been derived from the 
observed flux in the F160W filter since it samples the R-band in the 
rest-frame of the galaxies. To derive we used the relation 

M«(z) = F160W,„, - 5log[Di,(z)] - k^.^mwiz) (3) 

where F160W,o, is the total magnitude, Dl(z) is the luminosity dis- 
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tance [Mpc] at the redshift z of the galaxy and ^R,f leowCz) is the 
color k-correction term defined as 



kRJ^imW = \R{Z = 0) - F\6m{z)\,e,np 



(4) 



where the two magnitudes R{z = 0) and F 160^(2) arc derived from 
the best fitting template at z = and redshifted at the redshift of the 
galaxy respectively. The uncertainty afi'ecting this k-correction is 
typically comparable or smaller than the photometric errors since 
the filter F160W is extremely close to the rest-frame R-band over 
the whole redshift range considered and thus the dependence on the 
best fitting template tends to vanish. Indeed we have verified that 
even considering the oldest template among the best fitting tem- 
plates (the one 4.0 Gyr old) and the youngest template (1 Gyr old) 
the dilference between the k-correction is less than 0.15 mag over 
the whole redshift range considered. Thus, even hypothesizing to 
fail the fit to the observed SED of a galaxy we would wrong its ab- 
solute magnitude by no more than 0.1 mag. For this reason we con- 
sider our estimate of the R-band absolute magnitude of our galaxies 
extremely reliable. 

The stellar mass M, of the galaxies we derived is the one usu- 
ally computed in the literature and it is given by the equation 



M. = 2 • 10" 



• b ■ 4nDl(z) ■ MT'^'lLe 



(5) 



where b is the normalization factor of the best fitting model pro- 
vided by hyperz and >l™<"^«' is the mass associated to the best 
fitting template considering only the stellar mass still locked into 
stars. This quantity is listed in column 7 of the .4color files of CB08 
models. The mass we derived is the mass locked into stars at the 
epoch of their observation after the gas fraction returned to the in- 
terstellar medium. This mass is tj'pically about 60% the one derived 
including the gas return fraction, i.e. the one obtained by integrating 
the SFR over the age of the galaxy. A detailed comparison between 
different stellar mass estimators is given by Longhetti et al. (2008) 
which provide also the relations to convert an estimate to another 
accounting for difi'erent IMFs. The uncertainty affecting our stellar 
mass estimate depends mainly on the uncertainty affecting the SFH 
(t) and the age of the best-fitting model and on the best-fitting pa- 
rameter Ay, three parameters tightly linked among them. The SFH 
and the age affect mainly the value of M™"*"' while the extinction 
Ay afffect mainly the normalization b of the fit. Since we deal with 
galaxies of known spectral type and redshift the best-fitting SFH 
and age are sharply constrained producing negligible differences in 
the values of M™*'. For instance, even considering the youngest 
and the oldest ages possible in this range of redshift, i.e. 1 G3T and 
4 Gyr, the corresponding M™'"'"' for the same SFH would differ 
only by a factor 1.15. On the contrary, the normalization b of the 
model which is a free parameter in the fitting, can vary up to a fac- 
tor ~2 since it depends on the photometric accuracy in the various 
bands, on the number of photometric points sampling the spectrum 
of the galaxy and on the free parameter Ay. Thus, the internal ac- 
curacy of our stellar mass estimates is within a factor two and it is 
two times the uncertainty affecting the absolute magnitudes. This 
internal error does not consider the possible systematics due to dif- 
ferent IMFs (e.g. Salpeter IMF provides higher stellar masses than 
Kroupa and Chabrier IMFs) or different library models (see e.g. 
Maraston et al. 2006; Longhetti et al. 2008). Thus, for a compari- 
son with other samples such possible systematics should be taken 
into account by scaling, if necessary, the different estimates. Sys- 
tematics and scaling relations among different library codes can be 
found in Longhetti et al. (2008). 

The mean age of the best fitting model depends both on the 
SFHs and on the Ay in the way that higher Ay and shorter SFHs 
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Figure 3. The extinction Ay is plotted as a function of the mean age AgSgOl 
of the best fitting model for the sample of 32 ETGs. No systematic are 
present between the two parameters. 



provide younger best fitting models for a given observed spectral 
energy distribution. As previously said, most of the galaxies (24 out 
of 32) are best-fitted by the shorter SFHs considered, i.e. by models 
with T < 0.3 Gyr and Ay < 0.3 suggesting that the best fitting does 
not tend toward either old or young models, old with respect to the 
age of the universe at the redshift of the galaxy. However, the main 
degeneracy is between the age of the best-fitting model and the 
extinction Ay. Indeed, while two different SFHs, for instance the 
one with t ^ 0.1 and the other t < 0.3, accounts for differences in 
the best-fitting age of the order of few tenth of Gyr, different values 
of Ay even within the range <Ay < 0.6 can produce differences 
as large as 1 Gyr. An old stellar population can, in fact, be fitted 
by a young model reddened by an extinction Ay > 0. In order to 
verify the absence of any systematics in our results due to the best 
fitting procedure we have compared the values of the extinction Ay 
with the age of the best fitting template. In Fig. 3, we plot Ay as 
a function of the age of the best fitting model for the sample of 
32 ETGs. No systematics are present among the two parameters 
confirming that yoimg best fitting templates are not a faked result 
of the fitting procedure. It is worth noting that, with the exception of 
two galaxies of the GDDS sample, our estimate of the ages agrees 
within ~ 0.3 Gyr (0.5 Gyr for the old ones) with the ages derived 
by the various authors from the spectral features of the galaxies. 

In Tab. 2 we report for each galaxy the R-band absolute mag- 
nitude, the stellar mass and the mean age derived for each galaxy 
from the best fitting to the photometry. 



4 THE EVOLUTION OF THE KORMENDY RELATION 

The Kormendy relation (KR, Kormendy 1977) is a linear scaling 
relation between the logarithm of the effective radius R^ [Kpc], i.e. 
the radius containing half of the Ught, and the mean surface bright- 
ness (ju)e [mag/arcsec^]: 
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{juh = a+/3log(,R,) (6) 

The ETGs follow this relation with a fixed slope /3 ~ 3 up to 
z ~ I (e.g. Di Serego et al. 2005) while the zero point a varies with 
the redshift reflecting the evolution that the galaxy underwent. 

In Fig. 4 (left panel) the values of and </j>f in the R band 
for the 32 galaxies at 1 < z < 2 derived as described in the previous 
section are plotted on the |/<<,>Re] plane. The values of <yu)f have 
been corrected for the cosmological dimming factor (1 + z^. We 
converted the surface brightness in the F160W filter into 

that in the rest-frame R-band (/i>f applying to each galaxy the k- 
correction ^/(.fisovv described in the Sec. 3.2. In Fig. 4, the observed 
KR in the R band at z ~ (see La Barbera et al. 2003) 

18.2 + 2.92tos(^,), z = (7) 

is also shown (thin black line). The thick (red) solid line is the KR 
we obtain from the best fitting of our sample 

(ju), = 16.r«i + 2.nt°ilogiR,l l<z<2. (8) 

and the two thick dashed lines represent the dispersion at Icr of the 
relation. It can be seen that the slope of the KR we fit at z ~ 1.5 
does not change significantly with respect to the KR at z = while 
we detect a significant evolution of the zero-point a which changes 
more than 2 magnitudes in this redshift range. Thus, at least up to 
z ~ 1.5-2 ETGs tend to distribute on the [/Ue,Re] plane according 
to a KR with a slope similar to the one of local ETGs z = 0. The 
diflrerent zero-point accounts for the evolution which undergo the 
galaxies and tell us that in case of pure luminosity evolution, i.e. at 
constant R,,, galaxies must evolve by more than 2 magnitudes in the 
rest-frame R band from Zmean ~ 1.5 to z = 0, in agreement with pre- 
vious results (see e.g. Mcintosh et al. 2005; Longhetti et al. 2007; 
Cimatti et al. 2008). Such evolution exceeds almost 1 mag the one 
expected assuming an average passive luminosity evolution (PLE) 
for the whole sample and exceeds much more the evolution inferred 
from the observed luminosity function of galaxies in this redshift 
range (e.g. Feulner et al. 2007; Cirasuolo et al. 2007; Marchesini et 
al. 2007; Zucca et al. 2006; Pozzetti et al. 2003). However, our sam- 
ple of galaxies is spread over a range of redshift, I < z < 2, which 
spans about 2.5 Gyr of time, a large interval if compared to the cos- 
mic time at z ~ 1 (5.7 Gyr) and at z ~ 2 (3.2 Gyr). Thus, in order to 
assess whether the observed evolution of the zero point a of the KR 
can be accounted for by luminosity evolution, we derived for each 
galaxy its own luminosity evolution £(z) in the rest-frame R-band 
over the interval At = t{z = 0) - t(z), corresponding to the time 
elapsed from z to z = 0. The evolutionary term E(z) has been cal- 
culated taking into account the dilferent ages of the galaxies at the 
observed redshift as provided by the best-fitting model. Thus, for 
each galaxy we computed the term E{z) = [Mniz) - Mg(z = O)]mode/ 
i.e. the difference between the R-band absolute magnitude of the 
best fitting model with age Age and the R-band magnitude of the 
same model with age Age + At. This term added to equation (3) 
provides the surface brightness <jx}f that our galaxies would have 
at z = in the case of pure Ivmiinosity evolution. In Fig. 4 (right 
panel) we show how the 32 ETGs of our sample would be displaced 
at z = in the [/i(.,R(.] plane. We see that almost a half of the sam- 
ple reaches a surface brightness not exceeding the one derived from 
the KR at z = while the remaining half shows a surface bright- 
ness in excess of ~1 magnitude with respect to the local KR. In 
particular, the luminosity evolution brings 1 1 ETGs onto the KR 
at z = and two ETGs just below it while leaves the remaining 
19 ETGs with a surface brightness exceeding much more than one 
sigma the local KR. It seems that for a fraction of ETGs at z > 1 



the expected luminosity evolution is sufficient to account for their 
surface brightness. They move in the [jUf ,R(.] plane from high to low 
redshift in agreement with the local KR. On the contrary, for the re- 
maining fraction of ETGs the expected luminosity evolution is not 
sufficient to dim their surface brightness to the one defined by the 
local KR. The other parameter involved in the KR relation, the ef- 
fective radius R^, must evolve and the hypothesis of fixed size must 
be rejected for these galaxies. A size evolution of at least a factor 
~2.5 from z ~ 1.5 to z = is needed to account for the observed 
surface brightness excess. 



5 YOUNG V5 OLD ETGS: TWO DISTINCT 
POPULATIONS AT Z ~ 1.5 ? 

In the previous section we have seen that for 13 ETGs the pure lu- 
minosity evolution can move them from the KR at z - 1.5 onto 
the KR at z = 0, while for the remaining 19 ETGs a different more 
complex evolution is required. In fact, this suggests that two distinct 
populations of ETGs exist at z - 1.5. In order to better investigate 
this evidence, we compared the properties of the 13 ETGs whose 
luminosity evolution places them onto the KR at z = with the re- 
maining 19 ETGs. Basically, the luminosity evolution E(z) derived 
by models over a given interval At of time depends on the SFH and 
on the age of the best-fitting model. The SFHs which fit the SED of 
our galaxies are described by e-folding time t much shorter than the 
typical At over which the luminosity evolution is computed (t < 0.3 
Gyr to be compared with At ~ 9 Gyr from z = 1.5 to z = 0). Thus, 
the slightly different values of t cannot produce significant differ- 
ences in E{z) over this interval. Consequently, the other parameter 
affecting the luminosity evolution, that is the age of the best fitting 
model, must be the reason of the different behaviour of the two sub- 
samples. If so, we expect that the two sub-samples of galaxies show 
a different age distribution. 

In Fig. 5 the age distribution of the 13 galaxies (dashed green 
histogram) which agree with the KR at z = is compared with the 
age distribution of the remaining 19 ETGs (solid red histogram). It 
is worth noting that the two sub-samples, with the exception of 2 
ETGs (one for each sub-sample), describe two separated distribu- 
tions. The first distribution is sharply picked at ~ 1 Gyr, the other 
distribution picks at ~ 3.5Gyr. As shown in Sec. 3.2 age and extinc- 
tion Ay, are not correlated and no systematic are present between 
them (see Fig. 3). Moreover, the diff'erent values of t considered 
cannot account for such difl'erent ages. Indeed, a galaxy 1 Gyr old 
fitted by a SFH t = 0. 1 Gyr would be ~ 1 .3 Gyr old if fitted with a 
model T = 0.3 Gyr. Analogously, a galaxy 4 Gyr old fitted by a SFH 
T = 0.3 Gyr would be ~ 3.5 Gyr old if fitted with a model t = 0.1 
Gyr. Thus, the two different distributions are not a consequence of 
the degeneracy between SFH, age and extinction but they reflect 
real differences among the ETGs: the 13 ETGs which fall on the 
local KR are, in fact, younger than the remaining 19 ETGs. Given 
the uncertainties discussed above, it is reliable to consider a mean 
difference of about 1.5-2 Gyr between the age of the two popula- 
tions. Hereafter, we will refer to these two populations as young 
ETGs (yETGs) and old ETGs (oETGs). 

In fact, the different age of the stellar populations of the two 
sub-samples of ETGs is the reason of their different behaviour with 
respect to the KR. Figure 6 displays the evolution term E{z) in the 
R-band as a function of the cosmic time (bottom x-axis) or redshift 
(upper X-axis) for two galaxies whose stellar populations at z = 1.5 
are 1 Gyr old (dashed red Une) and 4 Gyr old (soUd red Une). The 
thin lines refer to a SFH with t = 0. 1 Gyr while the thick lines refer 
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Re fkpcl [kpc] 

Figure 4. Left panel - Mean surface brightness versus effective radius R,, [kpc] for the 32 ETGs of the sample. The thin (black) solid line represents the 
KR at z ~ (eq. 7) and the short-dashed lines represent the ±lo" dispersion of the relation. The thick (red) line is the KR at z ~ 1.5 as resulting from the best 
fit to the sample (eq. 8). The short dashed lines represent the scatter around this relation. All the data have been corrected for the cosmological dimming factor 
(1 + ?)■* thus, the deviation from the KR at z = reflects the evolution of the SB due to the luminosity and/or size evolution of galaxies. Black points are the 
S2# ETGs (Longhetti et al. 2007), cyan squares are the SA# ETGs from the GDDS (McCarthy et al. 2004), blue pentagons are the CIG# ETGs in the RDCS 
0848+4453 (Stanford et al. 1997), green hexagon are the HDF# ETGs, red triangles are the HUDF# ETGs (Cimatti et al. 2008) and the red starred symbol is 
53W091 (Dunlop et al. 1996). Right panel - Mean surface brightness versus effective radius R, [kpc] for the 32 galaxies of our sample after than they 
have evolved in luminosity. Each galaxy has evolved in luminosity over the interval At corresponding to the time elapsed from the redshift z. of the galaxy to 
Z = 0. The evolution E{z) has been derived accounting for the different age of each galaxy at the observed redshift (see Sec. 4). 
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Figure 5. Age distributions of the two sub-sample of ETGs. The dashed 
(green) histogram represents the distribution of the 13 ETGs whose surface 
brightness agrees with the local KR in case of luminosity evolution. The 
solid (red) histogram represents the distribution of the remaining 19 ETGs 
whose SB exceeds by more than one <t the KR. 
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Figure 6. Evolution term £(z) as a function of the cosmic time (bottom x- 
axis) or redshift (upper x-axis) for two galaxies which at z = 1.5 are 1 Gyr 
old (dashed green line) and 4 Gyr old (solid red lines) respectively. The thin 
lines have been obtained with a SFH described by t = 0.1 Gyr, the thick 
lines refer to a model with t = 0.3 Gyr. 
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to a r = 0.3 Gyr. In fact, we see that the difference in the luminosity 
evolution E{z) between the old and the young stellar population is 
about 1 mag at z = independently of the different values of t. 
Thus, the different mean age of the stellar populations of ETGs at 
z ~ 1.5 is the reason of their different expected evolution. We have 
compared also the absolute magnitude distribution and the stellar 
mass distribution of these two poptilations to gather other informa- 
tion about their evolutionary status and to search for signs of differ- 
ent history of star formation and mass assembly. In Fig. 7 we show 
the distributions of the absolute magnitude (left panel) and of the 
stellar mass (right panel) of the yETGs (dashed red histogram) and 
of the oETGs (solid red histogram). We see that yETGs tend to be 
less luminous and, accordingly, less massive than the oETGs even if 
the effect is not statistically significant, as confirmed by the KS-test 
performed to compare the two distributions (P(D > D„ax) = 0.02). 
However, it can be seen that the high-luminosity/mass tail is popu- 
lated only by old ETGs while the low-luminosity/mass tail is com- 
posed of only young ETGs. 

The SED fitting of our ETGs is based on optical and near-IR 
photometry which at z ~ 1.5 samples the wavelength range A < 0.8 
fim. In practice, at this redshift, 6 out of the 8-9 photometric points 
sample the UV and blue rest-frame emission of the galaxies whose 
continuum shape is affected by star formation episodes even if in- 
volving a negligible fraction of the stellar mass (see e.g. Fig. 2 in 
Longhetti et al. 2008). For this reason, it is more appropriated to 
consider the age we derived as a lower limit to the time elapsed 
since the last episode of star formation. If this latter is the major 
one, then this age will provide the formation redshift of the stel- 
lar population. Given the different age of the two populations of 
ETGs it is likely that the formation redshift of their stellar popula- 
tions are different. In particular, the stellar component of the oETGs 
{Age„„d ~ 3.5) formed at z/ ~ 5 - 6 while the stellar populations of 
yETGs {Age„,gci ~ 1) formed at z/ > 2.5 - 3, or at least the youngest 
population. 



6 THE SIZE-LUMINOSITY AND THE SIZE-MASS 
RELATIONS 

The dependence of the characteristic size of galaxies on their lumi- 
nosity and on their stellar mass has been recently studied for large 
samples of local galaxies (e.g. Shen et al. 2003), at intermediate 
redshift (z ~ 1, Mcintosh et al. 2005) and at high-redshift (z < 3) 
(e.g. Trujillo et al. 2004, 2006b; Cimatti et al. 2008; van Dokkum 
et al. 2008). Shen et al. (2003) on the basis of the SDSS data, show 
that the size of local early types and late types increases, as ex- 
pected, according to their luminosity and stellar mass. However, 
the ETGs follow steeper relations than late types showing that the 
size of a stellar system is not simply a function of its mass and that 
the history of its mass assembly can affect these relations (Shen et 
al.2003). Mcintosh et al. (2005) study the evolution of these rela- 
tions up to z < 0.8 - 1 by combining GEMS data (Galaxy Evolution 
from Morphology and SEDs, Rix et al. 2004) with COMBO-17 
data (Classifying Objects by Medium-Band Observations, Wolf et 
al. 2003, 2006b) while Trujillo et al. (2007) extend the study up to 
z ~ 2 thanks to the DEEP-2 survey (Davis et al. 2003; Bundy et al. 
2006) and to z ~ 3 (Trujillo et al. 2004, 2005) thanks to the FIRES 
data (Franx et al. 2000). Our data allow us to define these relations 
for a sample of ETGs at 1 < z < 2 with secure spectroscopic red- 
shift and classification and on the basis of the morphology derived 
by their red rest-frame continuum. 

The size-luminosity (S-L) relation for our sample at z ~ 1.5 



is shown in Fig. 8 (left panel) and it is compared with the relation 
found in the r band by Shen et al. (2003) using the Sersic half light 
radius for local ETGs (solid Une) 

log(R,) = -0.26Mr - 5.06 (9) 

The dotted lines represent the scatter of the relation. Young 
ETGs are marked with (green) squares while old ETGs with (red) 
points. The R-band absolute magnitude Mr of our galaxies is the 
one at the redshift of the galaxies. The offset with respect to the 
relation of Shen et al. reflects the evolution which ETGs undergo 
from their redshift to z = 0, the same evolution observed in the 
comparison of the KR shown in Fig. 4 (the one at z ~ 1.5 and that 
at z = 0). It is worth noting that the difference between the r pho- 
tometric band of the SDSS and the R Cousins band we use is about 
0.2 mag (Fukugita et al. 1995). Thus the possible uncertainties re- 
lated to the transformation between the two filters are negligible. 

In the right panel of Fig. 8 we show how the 32 ETGs of our 
sample would be displaced at z = in the [Mr, R^] plane in case 
of pure luminosity evolution, i.e. the R-band absolute magnitude 
plotted is Mr(z = 0) = MrIz) - E(z). In this case the different 
behaviour shown by the young ETGs with respect to the old ETGs 
is even sharper than in the case of the KR (Fig. 4). It is evident the 
agreement between the young ETGs and the local size-luminosity 
relation once considered their own luminosity evolution. All the 
yETGs are located within the scatter region of the z ~ relation. On 
the contrary, it is evident the disagreement between the old ETGs 
and the local S-L relation, disagreement which shows clearly that 
R(, must change from z ~ 1.5 to z = 0. It is worth noting that the 
calibration of the S-L relation based only on the oETGs data is: 

log(R,) = -0.28Mk(z = 0) - 5.6 

(dashed (red) line in Fig. 8), a relation with the same slop of the one 
at z ~ but with an offset of ~ 0.5 in the zero point corresponding 
to a factor ~ 3 in R^. 

In Fig. 9 our galaxies (filled symbols) are plotted on the size- 
mass (S-M) plane and compared with the S-M relation found by 
Shen et al. (2003) for the local ETGs (solid line) expressed by the 
following equation 

Re = 2.88 X 10-\M,/Mef^'' (10) 

The stellar mass they use is the one from Kauffmann et al. (2003) 
based on the Bruzual and Chariot (2003, BC03) models and on the 
Kroupa (2001) initial mass function while we used the CB08 mod- 
els and Chabrier IMF. Longhetti et al. (2008) show that the stellar 
mass obtained with Kroupa IMF differs by less than 5% from the 
one obtained with Chabrier IMF and that the use of BC03 mod- 
els leads to over-predict the mass estimate by a factor I.2-I.3 with 
respect to CB08 models. Thus, we decided not to apply any scal- 
ing factor to the relation found by Shen et al. (2003) given such 
small diff erences. Figure 9 shows that while 9 out of the 13 yETGs 
(70%) follow this size-mass relation at z ~ only 4 out of the 
19 oETGs (20%) agree with this relation. The yETGs for which 
the luminosity evolution is requested to bring them on the local 
S-L relation, do not need any mass or size evolution since they 
naturally match the local S-M relation. In other words, this is an 
evidence that yETGs are fully compatible with a simple evolution 
from z ~ 1.5 - 2 to z = of their At/L ratio due to pure luminosity 
evolution while their stellar mass remains unchanged. On the con- 
trary, the old ETGs follow well defined S-L and S-M relations but 
almost all of them have sizes much smaller than at z ~ implying 
that they have changed significantly their structure from z ~ 1.5-2 
to z = 0. Using the relation of Shen et al. (e.q. 10) we have derived 
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Figure 8. Left panel - Size-luminosity relation for our sample of yETGs (green squares) and oETGs (red points) at z ~ 1.5 compared with the local relation 
found by Shen et al. (2003, solid line). The dotted lines represent the scatter of the relation. Right panel - Size-luminosity relation in the case of pure luminosity 
evolution. It is shown how the 32 ETGs of our sample would be displaced at = in the [R, .M/;] plane in case of pure luminosity evolution, i.e. the R-band 
absolute magnitude plotted is Mr(z = 0) = Mr(z) - E(z). The dashed (red) line is the best fitting relation logR^ = -0.28M«(z = 0) - 5.6 to the oETGs (see 
Sec. 6). 



the mean value = l/A'Z(^f.o/^f,c), i-e. the ratio between the 
radius R,, of the 15 oETGs which do not follow the local S-M re- 
lation and the radius of the local ETGs with similar mass. We have 
obtained (/r^) = 2.6 ± 0.5 in agreement with the value 3.4+1.7 
found by Cimatti et al. (2008) in the same redshift range. Thus, 
the oETGs must increase their size by a factor 2.6 from z ~ 1.5 to 



z = 0, consistently with the result obtained from the KR and the S-L 
relation. It is worth noting that even hypothesizing that we have sys- 
tematically overestimated by a factor 2 the mass of all the oETGs 
their effective radius would be still 1.7 times larger then locally, i.e. 
they would be 5 times denser. To move the oETGs onto the local 
S-M relation their stellar mass should be 6 times smaller, a factor 
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Table 2. Morphological parameters of galaxies. The effective radii obtained by fitting the Sersic profile and the resulting surface brightnesses take into account 
the correction for the underestimate of derived from the simulations (see Sec.3. 1). We appHed a correction of 0.07 arcsec to the NIC2 data and of 0.03 arcsec 
to the NIC3 data. The error on the effective radii takes into account both the formal error of the profile fitting and the rms observed in recovering the intrinsic 
radius of simulated galaxies, namely a^'^ = 0.02 arcsec and crjy"'^ = 0.04 arcsec. The magnitude F160W(o( is the total magnitude of the best fitting profile 
as derived by galfit. The effective radius r^*^" is the original effective radius estimated by other groups. 
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too large to be accounted for by any model assumption. In Fig. 9 
a collection of ETGs at 1.3 < z < 2 (open symbols) taken from 
the literature is also shown. Triangles, circles and squares mark the 
GMASS galaxies from Cimatti et al. (2008), the two HUDF galax- 
ies from Daddi et al. (2005) and the two spheroids of McGrath et 
al. (2008) respectively. The redshift and the spectral type of these 
ETGs are spectroscopically confirmed and the morphology is based 
on HST-ACS observations. According to the analysis performed by 
the authors and to the parameters they derived, we divided this sam- 
ple of ETGs in old and young defining old those ETGs with age 
larger than 2 Gyr. The old ETGs are marked by red open symbols 
while the young ETGs are marked by green open symbols. It is re- 
markable the agreement with the behaviour shown by our sample 
of 32 galaxies: yotmg ETGs tend to distribute according to the lo- 
cal S-M relation (6 out of 9) while none of the old ETGs follow the 
local S-M relation. 

It should be noted however that a non negligible fraction of 
yETGs, both in our sample of 32 ETGs and in the sample taken 
from the literature, does not follow the S-M relation but follows the 
relation defined by the old ETGs. Thus, for some yETGs a major 



size evolution is still required. If this result will be confirmed on a 
more solid statistical ground proving that it is not due to an internal 
scatter in the estimate of the physical parameters (age, stellar mass 
and effective radius), it implies that yETGs follow different histo- 
ries of assembly and are less homogeneous than old ETGs. We will 
try to constrain their evolutionary path in the next section. 



7 CONSTRAINING THE FORMATION AND THE 
EVOLUTION OF ETGS 

The analysis performed in the previous sections shows that two 
populations of ETGs exist at z ~ 1.0- 2.0. They differ substan- 
tially for the age of their stellar populations by about 2 Gyr and for 
the scaling relations they follow. It is natural to ask how these two 
populations evolved from z < 2 to z ~ to match the properties of 
the local ETGs and which is their assembly history they followed 
to have the properties shown at z ~ 1.5. We have tried to answer 
to these questions placing our results in the hierarchical paradigm 
of galaxy formation and evolution taking into account the results 
obtained from various renditions of merging models. 
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Figure 9. Size-mass relation for our sample of ETGs (filled symbols) com- 
pared with the local relation found by Shen et al. (2003, solid line). The dot- 
ted lines represent the scatter of the relation. Red symbols mark the oETGs 
while green symbols mark the yETGs The two vectors represent the relation 
found by Boylan-Kolchin et al. (2006) for the two extreme values a = 0.6 
and a = 1.3. The open symbols (triangles, circles, squares) indicate the 
galaxies of the GMASS sample (Cimatti et al. 2008), the two passive galax- 
ies in the HUDF from Daddi et al. (2005) and the two spheroids of McGrath 
et al. (2008). We classified as old (and consequently marked with red color) 
those galaxies older than 2.0 Gyr. 



7.1 Tracing the evolution at z < 2 

The older ETGs of our sample at z 1.5 - 2 do not follow the S-M 
relation of local ETGs as well as the other scaling relations. Pure 
luminosity evolution from their redshift to z = does not bring 
them onto the local KR and S-L relation. oETGs are characterized 
by effective radii ~2.5-3 times smaller than those of the local 
ETGs with comparable surface brightness, absolute magnitude and 
stellar mass as deduced from the comparison with the local KR, 
S-L and S-M relations. Thus, an evolution of their size between 
z ~ 1.5-2 and z = must occur to bring them onto the local 
scaling relations. Such size evolution is often used to advocate the 
merging processes the ETGs should experience during their life in 
the hierarchical paradigm of galaxy formation and evolution (e.g. 
Trujillo et al. 2004, 2007; Bell et al. 2006; De Lucia et al. 2006; van 
Dokkum et al. 2008; van der Wei et al. 2008). Merging is indeed 
usually invoked as the most obvious and efficient mechanism to 
increase the size of galaxies. 

Boylan-Kolchin et al. (2006, BK06 hereafter), using simula- 
tions of dissipation-less merging, the so called "dry merging" (e.g. 
van Dokkum et al. 2005; Bell et al. 2006 and references therein) 
show that the remnants of dry mergers lie on the fundamental plane 
(FP, Djorgovsky and Davis 1987; Dressier et al. 1987) of their pro- 
genitors. However, the locations of the remnants in the projections 
of the FP, in particular on the R^-AI, relation, depends strongly on 
the merger orbit. Thus, the projections of the FP can provide a tool 
to investigate the assembly history of ETGs. In their analysis, they 
find that the expected increase of the size of an ETG due to merging 
follows the relation R^ oc M", with 0.6 < a < 1.3 (represented by 



vectors in Figure 9) depending on the orbital properties (see also 
Nipoti et al. 2002; Ciotti et al. 2007). They show also that the in- 
dex a is almost independent of the mass ratio of the progenitors 
(see also Khochfar and Silk 2006 for a similar result) suggesting 
that their findings are applicable both to minor and major mergers. 
We have tried to consider this model of dry merging to increase the 
size of oETGs. We have seen that the effective radii of oETGs must 
increase by a factor (fR^) — 2.6 from z ~ 1.5 to z = in order 
to match the local S-M relation. Thus, the condition Rf ^ 2.6i?, 
where Rj and Rf are the radii before (initial) and after (final) the 
merging must be satisfied. From the relation of BK06 it follows 
that M" - 2.6M" where Al, and Mf are the masses before and af- 
ter the merging. Consequently the mass Mf that the remnant must 
reach to increase the size 2.6 times is 



Mf = i.e^'^M, 



(11) 



The most efficient way to move oETGs from their location onto 
the local S-M relation is for o- = 1.3 (see vectors in Fig. 9), the 
maximum value found by BK06 which, by the way, minimizes the 
stellar mass of the remnant. We thus obtain 



Mf>2.lMi^ 



£»■< 1.3 



(12) 



i.e. the mass of the rerrmant is at least twice the mass before the 
merging. Any value of a lower than 1.3 would produce larger 
masses. This result is difficult to reconcile with the number density 
of high-mass ETGs in the local universe. Indeed, this mechanism 
would produce too much ETGs with masses much larger than 10" 
Mq and an evolution in the stellar mass density at z < 2 which is not 
observed (see e.g. Conselice et al. 2007). For instance, Saracco et 
al. (2005) show that the 7 galaxies S2F1# also studied in the present 
paper, account for 70% of the local population of ETGs with com- 
parable luminosity and mass. If they twice their mass/luminosity, 
at z ~ we should observe 2-3 times more ETGs with masses 
M, > 4 - 5 X 10" Mq then those in fact observed. Moreover, as 
previously noticed by Cimatti et al. (2008), it is difficult to imagine 
that given all possible orbital parameters in merging events, the ef- 
fective value of a is always close to the maximum one. Finally, we 
recall that values a < 1.3 would worsen the disagreement with the 
local number of high-mass ETGs. Thus, we conclude that merging 
cannot be the mechanism with which oETGs increase their size at 
z < 2 and that it is not the way to solve the problem. Other mech- 
anisms able to increase the size but to leave nearly unchanged the 
mass of ETGs must occur. Close encounters or, more generally, in- 
teractions between galaxies can act in this way. Their frequency 
and thus their efficiency depend on the number of close encoun- 
ters that a galaxy can experience in the last 9-10 Gyr of its life, a 
number that perhaps can be constrained from the statistics of pairs 
and from simulations. Minor or "satellite" merging (e.g. Naab et al. 
2007), i.e. merging between galaxies with masses Mi and M2 in 
the ratio ~0.1;1 or lower, would produce remnants with masses of 
the order of M2 but with larger size. It is not clear the ability of 
this kind of merging in enlarging the size, however it could act in 
the right way contributing to solve the problem of the small sizes 
of oETGs. 

The younger population of ETGs, the yETGs, follows the size- 
mass relation of local ETGs, with few exceptions. Luminosity evo- 
lution from z ~ 1.5-2 to z = 0, i.e. for fixed size R,, would bring the 
yETGs onto the local KR and S-L relation. Thus, for these galaxies, 
the evolution of the M/L ratio due to the expected luminosity evo- 
lution explains their observed properties at z ~ 1.5-2 and brings 
them to agree with the scaling relations of local ETGs. This sug- 
gests that the build-up of yETGs was already completed at z ~ 2 
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providing no evidence in favor of merging at z < 2 since it would 
bring them out of the S-M relation and the other scaling relations. 
Indeed merging, to move yETGs along the S-M relation, should 
take place for values a ^ 0.6 (see vectors in Fig. 9) producing 
remnants with masses Mf ~ 5At,. For the analogous reasons dis- 
cussed above, it is difficult to imagine that given all possible orbital 
parameters in merging events, the effective value of a for yETGs 
is always close to the minimum one. Moreover, the luminosity at 
z ~ of the remnant should dim according to the increased size in 
order to match the KR and the S-L relation. These requirements and 
fine tuning make this picture rather unlikely. Finally, the reasoning 
relevant to the exceeding number of remnant ETGs with masses 
well in excess to 10" M© applies also in this case, and corrobo- 
rates the conclusion that the assembly of yETGs was completed 
at z ~ 2 and that no merging has happened at z < 2. We cannot 
rule out that yETGs may experience satellite merging at z < 2 if it 
leaves nearly vmchanged their size besides their stellar mass. 

For the reasons here discussed we can conclude that ETGs, 
both young and old, have already reached their final stellar mass at 
z < 2. Major merging at redshift z < 2, if any, must necessarily in- 
volve a negligible fraction of the old ETGs while satellite merging 
could involve both yETGs and oETGs even if in a difi'erent way. 

7.2 Constraining tlie patli at z > 2: toward tlie formation of 

ETGs 

The older ETGs are characterized by a median age of about 3.5 
Gyr (and a dispersion of about 1 Gyr) which implies that their stars 
formed at z/ ~ 5-6. Given the short time they have at disposal 
to form masses of the order of 10" Mq of stars the star forma- 
tion rate was necessarily S FR > 100 Mg/yr (see also Cimatti et al. 
2008, 2004; McCarthy et al. 2004; Daddi et al. 2005; Longhetti et 
al. 2005; Feulner et al. 2005; Kriek et al. 2006). Moreover, oETGs 
are 2.5-3 times smaller than those at z ~ thus the physical mech- 
anism(s) acting at z > 2 must be capable also to produce very com- 
pact galaxies with stellar densities 15-30 times higher than the lo- 
cal ones and than the yETGs. Dissipational gas-rich merging can 
produce highly compact massive early-type galaxies if a high frac- 
tion of stars of the renmant formed during the merger in a violent 
starburst (e.g. Kochfar et al. 2006; Springel and Hernquist 2005). 
However, it is not clear whether the typical time scale of major 
merging can fit with the requirements above {Tmerge > 3 Gyr from 
Boylan-Kolchin et al. 2008). Naab et al. (2007), with their ab initio 
hydrodjmamics simulations, show that the early formation phase of 
galaxies can start with an initial burst of star formation at z ~ 5 
accompanied by mergers of gas-rich small sub-components and in 
situ intense star formation. Although some fraction of the stars is 
accreted (the "quiescent" component, Khochfar et al. 2006b) Naab 
et al. show that this phase has the characteristics of a dissipative 
collapse since it happens on very short timescale. Thus, the com- 
pactness of ETGs produced through the gas-rich merging of Koch- 
far et al. (2006) and the rapid dissipative collapse of Naab et al. 
(2007) provide a scenario which seems to fit, at least qualitatively, 
the intense SFR required at z ~ 5 - 6 and the need to assemble 
ETGs which at z < 2 are compact and old. 

The younger population, the yETGs, is characterized by a me- 
dian age of ~1 Gyr and a dispersion of about 0.2 Gyr which push 
the last burst of star formation at z/ > 2.5. In this case the con- 
straints on the possible physical mechanism(s) acting at z > 2 are 
less stringent than those required to assemble oETGs at the same 
redshift. A scenario in which yETGs accreted stellar mass by sub- 
sequent episodes of merging and by star formation in situ (which 



can also be triggered by satellite and secondary merging events) can 
qualitatively fit the properties of these yETGs which at z < 2 must 
appear younger and already enlarged in their size if compared with 
the local ones. The few of them which appear much more compact 
than the others yETGs could be the results of the gas-rich merging 
scenario proposed by Kochfar et al. (2006) and Naab et al. (2007) 
characterized by a rapid dissipative collapse happened at z > 2.5. 

From the observational point of view, the two populations of 
ETGs should have different progenitors. oETGs must experience a 
phase of intense star formation at high-z (z> 5), yETGs can experi- 
ence this phase at lower redshift or experience subsequent episodes 
of star formation possibly triggered by satellite merging. In any 
case a different epoch of formation and assembly must character- 
ize the two populations of ETGs. This result qualitatively agrees 
with the model of Khochfar and Silk (2006) which indicates that 
the scatter in the size of sirrular present-day ellipticals is a result 
of their formation epoch, with smaller ellipticals formed earlier 
through mergers much richer of gas than the mergers assembling 
larger ellipticals. In Fig. 10 the predicted evolution of sizes for 
ETGs with respect to the sizes of their local counterparts from the 
model of Khochfar and Silk (2006) is compared to the observed 
ratio R^(z)/R^(z = 0) for our sample of 32 ETGs. This ratio has 
been obtained for each galaxy dividing the observed effective ra- 
dius by that derived by the S-M relation of Shen et al. (2003) for 
the same mass of the galaxy. As expected oETGs are preferentially 
located on the curve representing the largest ratio which is expected 
for the highest-mass ETGs since they should form much earlier in 
the model. On the contrary yETG arc preferentially located onto 
or even above the curves representing the minimum size evolution 
which is expected for the lower mass ETGs which should form 
later. It is worth noting that while this model reproduces very well 
the observed relation between size evolution and formation epoch 
of ETGs (the older the more compact/denser) the correlation with 
the stellar mass is less evident from our data and a larger sample 
would be needed to probe this issue. As to the progenitor candi- 
dates, as suggested by Cimatti et al. (2008), a possible population 
of progenitors could be the sub-mm selected galaxies seen at z > 3 
(e.g. Blain et al. 2002; Tacconi et al. 2006, 2008) whose character- 
istics could fit those of ETGs at 1 < z < 2. Even if it is difficult 
to identify ffie progenitors we can tray to constrain ffieir redshift. 
Given the properties of oETGs (old and compact at z < 2) and the 
typical time scale of merging {r„erge > 3 Gyr), we should expect to 
see oETGs till z ~ 3 - 3.5, the only difference should be ffie age of 
ffieir stellar population correspondingly younger. On ffie contrary, 
yETGs should appear quite different at z > 2.5, most probably in 
the phase of merging, or star forming and interacting wiffi offier 
galaxies. 



8 SUMMARY AND CONCLUSIONS 

We presented ffie morphological analysis of a sample of 32 ETGs 
at 1 < z < 2 with spectroscopic confirmation of their redshift and 
spectral type based on HST-NICMOS observations in the F160W 
filter. These 32 ETGs have been selected from different samples 
and surveys on the basis of their spectroscopic and morphological 
classification and are characterized by a multiwavelength coverage. 
The HST-NICMOS observations in the F160W filter have allowed 
us to derive the effective radius R^ and the mean surface brightness 
(ju)e of galaxies in the rest-frame R-band of the galaxies, less af- 
fected by morphological k-correction and star formation ffien ffie 
optical bands usually used in ffie previous works. Through ffie best 
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Figure 10. The ratio R(,(z)/R(,(z = 0) (filled symbols) obtained dividing 
the observed effective radius of our galaxies by that derived by the S-M 
relation at z = of Shen et al. (2003) for the same mass is compared to 
the predicted evolution of sizes (curves) from the model of Khochfar and 
Silk (2006). Symbols are as in Fig. 8. Dotted, solid, dashed and dot-dashed 
curves represent the predictions for the different ranges of stellar masses 
listed in the legend. 



fitting of their SEDs at known redshift we derived the R-band abso- 
lute magnitude, the stellar mass and the age for each of them. The 
main results of the analysis we performed can be summarized as 
follows. 

• The 32 ETGs of our sample at 1 < z < 2 are placed on the 

plane defining a relation with the same slope of the KR 
at z ~ but with a different zero-point which accounts for the evo- 
lution they undergo from z - 1.5 - 2 to z = 0. We do not see 
differences between the 6 ETGs in clusters and the other ETGs in 
agreement with other works (e.g. Rettura et al. 2008; Gobat et al. 
2008) even if the very low statistics we have do not allow us a de- 
tailed comparison of different environments. 

• The ETGs of our sample are composed of two distinct popula- 
tions which differ for the age of their stellar populations: the older 
population, the oETGs, has a median age of about 3.5 Gyr which 
implies that the bulk of stars formed at z/ ~ 5 - 6; the younger 
population, the yETGs, has a median age of about 1 Gyr and corre- 
spondingly Zf > 2.5 - 3. Even if the absolute values of these ages 
can be model dependent, the different age cannot be accounted for 
by any model assumption while it is reliable an age difference of 
about 1.5-2 Gyr for the two populations of ETGs. yETGs tend to be 
less luminous and correspondingly less massive than oETGs even 
if this tendency is more evident in the tails of the distributions. 

• yETGs follow the size-mass relation locally observed (Shen et 
al. 2003) and their expected luminosity evolution from z ~ 1.5-2 
to z = at fixed size R^ brings them onto the local KR and size- 
luminosity relation. Thus, their surface brightness and stellar mass 
density do not exceed those of local ETGs with comparable lumi- 
nosity and stellar mass, i.e. no size evolution is required at z < 2. 
These properties suggest that young ETGs have already completed 



the growth of their stellar mass at z ~ 2 being them onto the S-M 
relation and that they must evolve at z < 2 purely in luminosity to 
match the local KR and S-L relation. This provides no evidence in 
favor of major merging at z < 2 since it would bring them out of 
the local scaling relations. 

• oETGs do not follow the local size-mass relation since they 
have sizes 2.5-3 times smaller then those provided by the S-M re- 
lation at their stellar masses. Pure luminosity evolution from their 
redshift to z ~ is not sufficient to bring them onto the local KR 
and S-L relation. Also in this case the effective radii are 2.5-3 times 
smaller than the local ETGs with comparable surface brightness 
and absolute magnitude. Thus, an evolution of their size at z < 2 
must occur to reconcile the oETGs with the local population of 
ETGs. Major (dry) merging at z < 2 cannot solve the problem since 
it would produce too much ETGs with stellar masses > 5- 10" -10'' 
Mq and it should happens only under particular orbit conditions to 
move oETGs onto the local S-M relation. Other mechanisms able 
to increase the size and to keep constant the stellar mass of oETGs 
(e.g. satellite merging, close encounters) must be invoked. 

• The different properties shown by the yETGs and the oETGs 
at 1 < z < 2 imply different evolutionary paths from their formation 
to z < 2. oETGs are much more compact and hence denser than lo- 
cal ones (15-30 times denser than the local ETGs with comparable 
stellar mass and than yETGs) they are old with respect to the age 
of the universe at their redshift. Thus, their stellar mass must have 
formed at high-z (z ~ 5) following a sort of dissipative gas-rich col- 
lapse ables to form rapidly most of the stellar mass thus producing 
a compact old remnant at z < 2. This scenario is qualitatively fitted 
by the merging models of Khochfar et al. (2006) and Naab et al. 
(2007). Considering the typical time scale of merging {r^erge > 3 
Gyr; Boylan-Kolchin et al. 2008) and the age of the stellar popula- 
tion of oETGs, we conclude that oETGs must exist as they are till 
z > 3 - 3.5 with stellar populations correspondingly younger. Pro- 
genitors should be searched for among the population of galaxies at 
z > 3.5 - 4 as high-mass galaxies with intense star formation. Pos- 
sible candidates could be the sub-mm/selected galaxies (e.g. Blain 
et al. 2002; Tacconi et al. 2006, 2007) as suggested by Cimatti et 
al. (2008). 

The formation scenario for the yETGs seem to be rather different 
even if less constrained with respect to the one of the oETGs. The 
age of the yETGs implies that the last burst of star formation has 
taken place close to their redshift, i.e. z ~ 2.5. They are not denser 
than the local one and are placed on the local S-M relation. They 
have completed their stellar mass growth as the oETGs but in a way 
such that their size is larger. Major merging, satellite merging and 
close encounters coupled with star formation in situ can qualita- 
tively fit these requirements producing at z < 2 ETGs with a young 
component of the stellar population and sizes comparable to those 
of the local ETGs with similar stellar mass, surface brightness and 
luminosity. This population of yETGs should appears at z > 2.5 - 3 
as star forming and/or interacting galaxies. 



On the basis of the above results we believe that a key observational 
test would be the measure of the velocity dispersion of oETGs and 
yETGs since such quantity would unambiguously address the ques- 
tion whether the two populations are dynamically different thus 
providing unique constraints on the mechanism of their formation 
and on their size evolution at z < 2. 
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APPENDIX A: PHOTOMETRY OF THE SAMPLE 

In Table 1 the multiband photometry for the whole sample of early- 
type galaxies is reported. With the exception of the F160W-band 
magnitude which we have estimated on the HST-NICMOS images, 
the other magnitudes have been taken from the literature. 

The photometry of the two galaxies in the HUDF previously 
studied by Daddi et al. (2005) and by Cimatti et al. (2008) comes 
from the catalog of Coe et al. (2006). Namely, the B, V, R, z and 
J magnitudes are the magnitudes measured through the HST filters 
F435W, F606W, F775W , F850W and FUOW respectively. The 
Ks-band magnitude are taken from Cimatti et al. (2008) and comes 
from VLT-ISAAC observations in the Ks filter. 

For the sample of Longhetti et al. (2005), the magnitudes are 
those from the MUNICS survey whose optical filters are slightly 
different from the standard Kron-Cousins filters. A detailed de- 
scription of the MUNICS photometry and of the filter response is 
given by Drory et al. (2001). 

The seven-filters (B, V, R, I, z, H, K) photometry of the 10 
ETGs selected from the GDDS (Abraham et al. 2004) was orig- 
inally taken from the photometric catalogs of the Las Campanas 
Infrared Survey (LCIR survey McCarthy et al. 2001). The observa- 
tions of the fields, the filters used and the photometric information 
are described in Chen et al. (2002) and Firth et al. (2002). 

The photometry of the 6 galaxies at z ~ 1.27 belonging to the 
cluster RDCS 0848+4453 in the Linx Supercluster is taken from 
the paper of Stanford et al. (1997). The three galaxies CIG_135, 
CIG_108 and CIG_142 seem to coincide with the three galaxies 
studied by van Dokkum et al. (2003). The photometry is described 
in Stanford et al (1997) and subsequent observations are quoted in 
Mei et al. (2006). 

The photometry of the three galaxies in the HDF-N previously 
studied by Stanford et al. (2004) are partly taken from the HDF-N 
catalog (v2). The magnitudes in the filters B, V and I are in fact the 
magnitudes measured through the HST filters F450W, F606W and 
F814W. 

The photometry in the R, H and K bands of the galaxy 
53W091 are taken from the papers of Dunlop et al. (1996) and 
Spinrad at al. (1997). The I and I magnitudes are the magnitudes in 
the F814W and Fl lOW filters respectively taken from Waddington 
et al. (2002). 
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